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Abstract. Low-mass protostars are less luminous than expected. This luminosity problem is 
important because the observations appear to be inconsistent with some of the basic premises 
of star formation theory. Two possible solutions are that stars form slowly, which is supported 
by recent data, and/or that protostellar accretion is episodic; current data suggest that the 
latter accounts for less than half the missing luminosity. The solution to the luminosity problem 
bears directly on the fundamental problem of the time required to form a low-mass star. The 
protostellar mass and luminosity functions provide powerful tools both for addressing the lumi- 
nosity problem and for testing theories of star formation. Results are presented for the collapse 
of singular isothermal spheres, for the collapse of turbulent cores, and for competitive accretion. 



1. The Luminosity Problem 

Why don't protostars shine more brightly? In a seminal paper, Kenyon et al. (1990) 
identified this luminosity problem, developed an approach to treat it (the protostellar 
luminosity function), and proposed almost all the solutions to the problem that have 
subsequently been studied. The luminosity problem is simple to state: The accretion 
luminosity of a protostar is 
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where / acc is the fraction of the accretion energy that goes into radiation (Kenyon et al. 
took / acc = 1), m is the protostellar mass, and rh is the accretion rate. Stahler (1988) has 
calculated the protostellar radius for a variety of masses and accretion rates, including the 
thermostatic effects of deuterium burning; we find that r* ~ 2i? Q is a typical value for the 
harmonic mean radius. By comparing the number of embedded sources with the number 
of T Tauri stars, Kenyon et al. (1990) inferred a star formation time tf ~ (0.1 — 0.2) Myr 
in Taurus- Auriga; to build up a star of average mass (0.5 Mq) in this time requires an 
accretion rate m ~ (2.5 — 5) x 10 -6 M Q yr _1 . On average, the mass of a protostar will be 
about half its final stellar mass, which implies that the typical luminosity of a protostar 
that will become a star of average mass is L ~ (10 — 20)L Q . (The median luminosity 
of the average mass star is comparable to the median luminosity — see ^S] below.) The 
median luminosity in the Kenyon et al. sample is 1.6L Q . One statement of the luminosity 
problem is that the median protostellar luminosity is observed to be about an order of 
magnitude less than the expected value. 

Kenyon et al. (1990) also provided an alternative description of the luminosity problem: 
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Identifying the peak in the observed luminosity function of embedded sources at 0.3 L & 
as the luminosity of the lowest mass stars, which they took to be 0.1 Mq (keep in mind 
that this was prior to the discovery of brown dwarfs), and estimating the radius of 
these protostars as r* ~ 1 Rq, they inferred a mass accretion rate of only 10~ 7 M Q yr -1 . 
However, very general arguments (Stahler, Shu, & Taam 1980) indicate that the accretion 
rate due to gravitational collapse should be of order 
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There are solutions that give higher accretion rates than this, such as the Larson (1969)- 
Penston (1969) solution for the collapse of an isothermal sphere of constant density, but 
there are no solutions that give lower accretion rates, at least prior to the time that the 
accretion is affected by the outer boundary (Henriksen et al 1997). Since the observed 
temperature in molecular clouds is T ~ 10 K, this again leads to an order of magnitude 
discrepancy between observation and theory. 

Results from the recent Spitzer c2d survey of five nearby star-forming molecular clouds 
(Evans et al. 2009; the survey does not include Taurus- Auriga) confirm that protostars 
have low luminosities. Evans et al. (2009) classified the young stellar objects (YSOs) in 
the traditional class system, in which Class represents protostars with envelope masses 
greater than the protostellar mass; Class I represents embedded objects with masses 
greater than the envelope mass; Flat Spectrum objects represent a transitional class; 
and Class II corresponds to T Tauri stars. Dunham et al (2010) analyzed this sample 
and found an extinction-corrected median luminosity of 1.5 Lq. Including 350 /im data 
for many of the brighter sources, they obtained an extinction-corrected mean of 5.3 Lq. 

The luminosity problem is thus well established observationally. It is fundamental 
because observations appear difficult to reconcile with some of the the basic premises of 
star formation theory, that stars form via gravitational collapse (eq. II. 3[) and that they 
radiate the binding energy in the process (eg. 11.11) . As the discussion above illustrates, 
there are two related aspects to the problem: the observed luminosity appears to be less 
than that theoretically expected, and there is an excess of very low-luminosity sources. 



2. Proposed Solutions 

Kenyon et al. (1990) proposed a number of solutions to the luminosity problem, of 
which the two major ones are slow accretion and episodic accretion (discussed below). 
They also suggested that brown dwarfs could alleviate the luminosity problem by pro- 
viding sources with luminosities less than they considered. Since their work predated the 
discovery of brown dwarfs, they did not include them. It is now known that brown dwarfs 
constitute about 20% of stars (Andersen et al. 2008), and they do permit some sources 
to have lower luminosities. We include brown dwarfs in the models described below. 

There is an additional effect that reduces the luminosity problem: The hydromagnetic 
outflows observed from protostars extract kinetic energy from the accreting gas (McKee 
& Ostriker 2007), reducing the energy radiated. If half the energy lost by the disk is 
mechanical energy and this in turn is half the total potential energy, then / acc ~ | 
(Offner & McKee 2010). 
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2.1. Episodic Accretion 

Kenyon et al. (1990) suggested that the accretion onto the protostar (as opposed to infall 
onto a circumstellar disk) could be episodic, so that much of the protostellar mass would 
be accreted in short periods of time, with high luminosities. Such brief, high-luminosity 
accretion events could be associated with FU Ori outbursts, which have inferred accretion 
rates of m ~ 1CT 4 M Q yr^ 1 . Kenyon et al. (1990) pointed out that the 150 YSOs in 
Taurus Auriga that formed in the last 10 6 yr correspond to a total accretion rate of 
0.75 x 10~ 4 M Q yr -1 , comparable to that of a single FU Ori object; thus, a signficant 
fraction of the mass of a protostar might be acquired during FU Ori events. 

Hartmann & Kenyon (1996) refined this argument: They cited a rate of star formation 
within 1 kpc of the Sun of (5 — 10) x 10 -3 Mq yr" 1 , which is similar to the recent estimate 
of 8 x 10~ 3 Mq yr" 1 by Fuchs et al (2009). At that time there were 5-9 known FU Ori 
objects within 1 kpc; if each were to accrete at a rate of 1O _4 M yr -1 , then protostars 
could gain about 10% of their mass this way. They pointed out that this was a lower limit, 
since more FU Ori objects would most likely be discovered. That has indeed occurred: 
Greene et al. (2008) count a total of 22 known FU Ori objects, of which 18 are within 
1 kpc. Rounding off, we infer that 20 objects accreting at 10 _4 M Q yr -1 would account 
for 25% of the mass accreted by protostars within 1 kpc. 

To elaborate on this result, consider a simple model for episodic accretion, in which 
protostars are in a high-accretion state for a total time Aihigh and in a low-accretion 
state for the rest of the time, which is effectively the star formation time, Aii ow ~ tf. 
The fraction of the mass accreted during a high state is 

F high = mhi ^ h . (2.1) 

We shall adopt parameters appropriate for FU Ori outbursts, but in fact all episodes of 
high accretion could be included in F high . The total time spent in the high state is 

At high = ^p, (2.2) 

where Af P: high is the number of protostars in a high state in some volume of the Galaxy 
and TV* is the star formation rate there. For a mean stellar mass of 0.5 Mq, the local 
star formation rate cited above corresponds to a birthrate within 1 kpc of the Sun 
of 0.016 stars yr" 1 . If the number of FU Ori objects within 1 kpc is about 20, then 
Athigh = 1250 yr. Inserting Hartmann & Kenyon's estimate that mhigh — 10~ 4 M Q yr -1 , 
we recover the result given above, fhigh — 0.25. We also note that the duty cycle of FU 
Ori outbursts is small: AihighA/ ~ 10 3 /(10 5-6 ) < 0.01. This small value is consistent 
with the absence of any known FU Ori sources in the Evans et al. (2009) sample. 

There are two uncertain numbers that entered into this estimate of -Fhigh, the number of 
FU Ori objects, A/^highi and the accretion rate, rhhigh- Undoubtedly, more such objects 
will be discovered within 1 kpc of the Sun in the future; however, it should be noted 
that several of the bursting objects, such as L1551 IRS5, have luminosities, and therefore 
accretion rates, an order of magnitude less than the average (Hartmann & Kenyon 1996). 
Further study will also refine the observed value of the mean accretion rate. Hartmann 
& Kenyon (1996) infer an accretion rate rhhigh = 1-9 x 10~ 4 M Q yr -1 and a stellar radius 
r* = 5.9 Rq for FU Ori itself. However, after the rapid accretion ceases, the star will 
shrink back to its original size, releasing a comparable amount of energy. A mean value 
TOhigh — 1 x 10~ 4 M Q yr -1 thus seems reasonable in this case. 

Hartmann & Kenyon (1996) also noted that given 5 known outbursts in 60 years and 
a star formation rate of 0.01 stars yr -1 within 1 kpc of the Sun implies that there are 
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about 10 bursts per object. The more recent data cited above lead to a similar conclusion. 
Since Aihigh — 1250 yr, this means that a typical outburst lasts about 100 yr, consistent 
with the observationally inlerred liletime (Hartmann & Kenyon 1996). 

Finally, we note that a potential problem with the FU-outburst solution to the lumi- 
nosity problem is that these outbursts appear to be located preferentially in regions of 
low star- formation rates (Greene et al 2008). 

2.2. Slow Accretion (m < 10~ 6 Af© yr^ 1 ) 

As an alternative solution to the luminosity problem, Kenyon et al. (1990) suggested two 
ways to increase the ages of the protostars and therefore reduce the inferred accretion 
rate: (1) If the star formation rate decreased substantially over the lifetime of the T Tauri 
stars, then the observed number of Class I sources would require a greater lifetime; or (2) 
if the lifetime of the T Tauri stars were larger, then the protostellar lifetime would increase 
proportionately. They rejected the first possibility since it is difficult to understand how 
star formation can decelerate in 1 Myr when it is in a region with a crossing time of 
10 Myr. The second possibility is more plausible, since the estimated lifetime of T Tauri 
stars has increased in the intervening 20 years. In their recent analysis of protostellar 
lifetimes, Evans et al (2009) concluded that the lifetime of Class II sources is about 2 
Myr, which would imply a Class I lifetime of tf ~ (0.2 — 0.4) Myr in Taurus Auriga. 
The larger of these two values is favored by data in the c2d survey, for which Evans et 
al. (2009) found a Class I lifetime of 0.44 Myr. The total protostellar lifetime, tf, must 
also include the time spent in the Class phase, which they found to be ~ 0.1 Myr, so 
that the total inferred mean lifetime of protostars in these clouds is (tf) = 0.54 Myr. 
This long protostellar lifetime, together with the correction for non-radiative energy 
losses (/_acc = |) and the correction for unseen outbursts (also a factor of |) effectively 
resolves the luminosity problem: The accretion rate for a star of average mass is then 
m ~ 10~ 6 M© yr , which corresponds to a typical accretion luminosity from equation 
(II. lj) of about 2.2 L©, comparable to the observed median luminosity of 1.5 L©. 

What about the other aspect of the luminosity problem, which is that the inferred 
accretion rates are much less than expected theoretically? There are two physical ef- 
fects that reduce the accretion rate below that in standard models, namely, protostellar 
outflows and the finite size of the protostellar envelope. Bontemps et al. (1996) found 
that protostellar outflow rates vary as the envelope mass and showed that this could be 
understood if accretion rates declined exponentially with time: m = m exp(— t/t*) (sec 
also Myers et al. 1998), where the decay time is i» = m//rho and is the final proto- 
stellar mass (ie., the initial stellar mass). This exponential decline in the accretion rate 
does not capture the reduction due to the initial stage of protostellar outflows (Shibata 
k, Uchida 1985). It is thus plausible when these effects are included, current theories 
can be consistent with the inferred low accretion rates. However, this will be possible 
only if the high accretion rates associated with an initial Larson-Penston accretion phase 
(Henriksen et al. 1997, Schmeja & Klessen 2004) do not release a signficant amount of 
radiative energy. 

3. The Protostellar Luminosity Function 

3.1. Previous Work 

Any solution of the luminosity problem must explain the distribution of protostellar 
luminosities-the protostellar luminosity function (PLF). Kenyon et al. (1990) introduced 
the PLF and considered two cases: (1) the standard isothermal sphere (IS) case (Shu 
1977), in which the accretion rate is independent of mass so that the formation time, 
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tf, is proportional to the final mass, my; and (2) a case in which the accretion rate is 
proportional to m/, so that i/ is independent of m/. Fletcher and Stahler (1994a, b) 
extended this work in the IS case by determining both the protostellar mass function 
(PMF) and the PLF for pre-main-sequence stars. 

Dunham et al. (2010) carried out a detailed analysis of the joint distribution of proto- 
stellar luminosities and bolometric temperatures. Interestingly, they adopted an accretion 
rate of 4.6 x 10 -6 Mq yr _1 , much higher than implied by the protostellar lifetimes inferred 
by Evans et al. (2009). In one series of models, they included outflows, which entrained 
envelope material and resulted in a core star-formation efficiency of 0.3-0.5. The outflows 
were not intrinsically collimated, as in the study by Matzner & McKee (2000), so that 
outflow cavities eventually expanded to an opening angle of 90 deg, terminating accre- 
tion. The outflows led to substantial variations in the observed bolometric luminosity and 
temperature as functions of the inclination angle. Despite the reduction in the accretion 
rate, the mean luminosities were still too high. The best agreement with observation 
was obtained by assuming that most of the mass was accreted in FU Ori-type events 
with accretion rates of 1O _4 M0 yr _1 , although the resulting models spent much more 
time with bolometric temperatures above 10 3 K than is observed. These models required 
-Fiiigh ~ 0.8 (Dunham, private communication). This is consistent with observation only 
if the number of FU Ori sources within 1 kpc of the Sun were about 3 times higher than 
currently known (see eq. 12. 2|) . 



We follow McKee & Offner (2010) and Offner & McKee (2010) in describing the protostel- 
lar mass function (PMF) and the protostellar luminosity function (PLF). Let i/j(mf)d\nmf 
be the fraction of stars born in the mass range dmf, ip is thus the IMF. The PMF is 
the present-day mass function of the protostars, and it must be consistent with the IMF 
when the protostellar mass reaches its final value, m/. We denote the fraction of pro- 
tostars in the mass range dm that will have final masses in the mass range dm/ by 
ijjp2{m,mf)d\nmd\nmf, it is the IMF weighted by the time the protostar spends at a 
given mass, 



where (tf) is the average protostellar lifetime. We have assumed that the accretion rate 
is a function of m and m/, which can be readily generalized to allow for high and low 
accretion states. The PMF is the integral of this function over all possible values of the 
final mass, ranging from the lower bound on the IMF, mi, to the upper bound, m u , 
subject to the constraint m ^ mf\ 



It follows that the fraction of protostars in the mass range dm and the luminosity 
range dL is 



3.2. The PMF and the PLF 



ifj p2 (m,mf) 
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ty p 2(L, m) dhxmdhxL = ipp2(m, m/) dhimdhimf. 



(3.3) 



The PLF is then 




(3.4) 



(3.5) 
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where the limits of integration are such that mi ^ mj ^ m u and m rrif. The PMF and 
PLF depend on the history of the mass accretion rate, m(m, to/) (eq. I3.1[) and therefore 
on the theory of star formation. 



4. Testing Theories of Star Formation 

4.1. Accretion Histories 

We consider three different theories of star formation, plus one variant: 

* Isothermal Sphere (IS, Shu 1977) — Inside-out collapse of singular isothermal 
sphere: 

to — TOjg = 1.54 x 1CT 6 (T/10 K) 3 / 2 M Q yr" 1 . (4.1) 

The formation time is proportional to the mass, tf am, so this accretion model is valid 
only for low-mass stars (Shu, Adams & Lizano 1987). 



* Turbulent Core (TC, McKee & Tan 2002, 2003)— Inside-out collapse of a turbulent 
core: 

1/2 

/ ni 1 

TOTC 



m 

TO/ 



\ 1/2 

J to 3 / 4 , with totc oc S 3/4 . (4.2) 



The accretion rate increases with both the protostellar mass, to, and the final mass, to/. 
This model was developed for high-mass star formation, and is equivalent to the accretion 
rate in the Hennebelle & Chabrier (2008) theory of the IMF. 

* Two-Component Turbulent Core (2CTC) — Begins as isothermal accretion and 
evolves to turbulent accretion (similar to the TNT model of Myers & Fuller 1992): 

1/2 

(4.3) 
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* Competitive Accretion (CA, Zinnecker 1982, Bonnell et al. 1997) — Stars form in 
a common gas reservoir, usually accreting at the tidally limited Bondi rate, to oc to 2 / 3 , 
and all having the same formation time: 

2/3 

/// = toca 1 ) mf) with toca oc 1/ (free-fall time). (4-4) 



TO 
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4.2. Comparing the PMF and PLF with Observation 
We adopt a truncated Chabrier (2005) IMF: 

"log 2 (TO//0.2)" 



■(/'(to/) oc exp - 



2 x 0.55 2 



OC TO/ 



-1.35 



(to/ ^ 1 M Q ) 

(1 Mq < TO/ ^ TO„) 



(4.5) 
(4.6) 



The Evans et al. (2009) sample of Class II YSOs has about 400 objects. In this sample, 
9 YSOs are expected with masses exceeding 3M Q if m u ^> 3M Q (i.e., stars with masses 
much greater than 3M are possible). Since none are seen, we adopt an upper cutoff 
m u = 3A/ Q for the IMF. 

Figure rjja, illustrates the PMF (eq. 3.2) estimated with each of the accretion histories. 
As discussed by McKee & Offner (2010), the models predict very different mass dis- 
tributions. For example, in the isothermal sphere case more massive protostars spend a 
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m (Mo) L (U) 

Figure 1. Left: The PMF, ip p (m), for the four accretion models and the Chabrier IMF. 
Right: The PLF, ^ P (L), for the same models assuming no tapering of the accretion rate, with 
Fhigh = 0.25 and (if) = 0.56 Myr with the data from Dunham et al. (2010). Both panels adopt 
a cluster upper mass m u = 3 Mq. 



longer time accreting and consequently weight the distribution towards higher masses. In 
contrast, for competitive accretion, where all the protostars share the same protostellar 
lifetime, the significantly larger number of low-mass protostars shifts the PMF towards 
lower masses. Unfortunately, since we can't directly measure the protostellar masses, 
Figure QJi is not sufficient to observationally discriminate between the models. 

For a more direct comparison, Offner & McKee (2010) calculate the PLF using the 
predicted mass functions in combination with a stellar evolution model (see Offner et 
al. 2009 for details) . Figure [T]d shows the PLF for each accretion model plotted with the 
extinction corrected luminosities from Dunham et al. (2010). The model curves assume 
that 25% of the mass is accreted during unseen bursts (i*high = 0.25) and (tf) = 0.56 Myr. 
This is equivalent to applying both the slow accretion and variable accretion solutions 
and thus likely represents a lower bound on the predicted luminosities. For comparison 
with Dunham et al. 2010, we adopt an upper bolometric luminosity uncertainty of 50% 
and use the uncorrected bolometric luminosities to set a lower error bound. 

The distributions can be characterized by the mean, median, and standard deviation. 
For m u — 3 Mq and L mm = 0.05 Lq, the fiducial models have means in the range 2.5 
L (2CTC) - 3.6 Lq (CA), all a factor of -1.5-2 below the observed value: 5.3 ±f^ Lq. 
We find that the mean luminosities fall in a narrower range, 2.6 Lq (2CTC)- 3.4 Lq (IS), 
when the accretion rates taper off towards the end of the protostellar lifetime: 



rh = mo(m, rrif) 



1-|1 



(4.7) 



where mo(m,m/) is the untapered accretion rate. (Foster & Chevalier 1993 found that 
the accretion rate tapered off at late times in their ID calculations, and Myers et al. 1998 
included an exponential tapering in their models.) Only the non-tapered CA and tapered 
IS models fall within the uncertainty, suggesting that the adopted lifetime may be too 
high by as much as a factor of 2. In contrast, the fiducial medians, which range from 0.8 
Lq (CA) to 2.5 Lq (IS), are in better agreement with the observed median of 1.5 i§=| Lq. 
The median of the TC (0.9 Lq) and 2CTC (1.4 Lq) models are within error and remain 
so even for a lifetime reduced by a factor up to 2.4 and 1.6, respectively. The observed 
standard deviation of log L, 0.7^-f , is consistent with the 2CTC (0.6), TC (0.7) and CA 
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(0.8) cases. The outcome of the comparison is sensitive to the values of Fhigh and / acc , in 
addition to (tf), and all have significant uncertainty. The parameter dependence may be 
reduced by comparing to the ratio of the median to the mean luminosity, 0.3 z^j- This 
rules out the IS model. 

5. Conclusions 

The luminosity problem in low-mass star formation can be resolved if low-mass stars 
form slowly, over a period ~ 0.5 Myr, as suggested by the results of Evans et al. (2009). 
Indeed, if some of the accretion energy is released mechanically (1 — f acc — 1/4) and in 
unseen FU Ori outbursts (-Fhigh — 1/4), then for most models the star formation time 
must be somewhat less than 0.5 Myr to be consistent with the observed protostellar lumi- 
nosities. Different theories of star formation predict different prototellar mass functions, 
which are currently inaccessible to direct observation, and protostellar luminosity func- 
tions, which have been observed. The latter serves as an important metric to discriminate 
between theories of star formation. 
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